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Abstract 

This note is concerned with the study of the initial boundary value problem for systems 
of conservation laws from the point of view of control theory, where the initial data is 
fixed and the boundary data are regarded as control functions. We first consider the 
problem of controllability at a fixed time for genuinely nonlinear Temple class systems, 
and present a description of the set of attainable configurations of the corresponding 
solutions in terms of suitable Oleinik-type estimates. We next present a result concerning 
the asymptotic stabilization near a constant state for general n x n systems. Finally 
we show with an example that in general one cannot achieve exact controllability to a 
constant state in finite time. 
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1 Introduction 



Consider an n x n system of conservation laws on a bounded interval 
u t + f(u) x = i>0, xe}a,b[, 

with the initial condition 

u(0,x) = <p(x), a < x <b, 
and a weak form of the Dirichlet boundary conditions 

Ui(t, a) = cti(t), Ui(t,b) = (3i{t), t>0 



(1.1) 
(1.2) 

(1.3) 



(see |Ol Q |l9| and reference therein for several weak formulations of (1.3)). 

We want to study the effect of boundary conditions on the solution of (IT) from the 
point of view of control theory. Namely, following the same approach adopted in Q for 
scalar conservation laws, we take the initial data tp fixed, and, regarding the measurable maps 
on, Pi in (1.3) as control functions, we want to investigate the property of the attainable set 
for (§"]])- (|T|), i.e. of the set 

A(T) = |u(T, •) ; u is a sol. to (1.1) - ([T|)} , (1.4) 

which consists of all profiles that can be attained at a fixed time T > 0, by entropy weak 
solutions of (1T)-(L3). For the definitions and the basic properties of weak solutions we 
refer to §. See also Jf| || 0, § for results concerning the existence and L 1 stability of 
entropy weak solutions of the mixed problem taking values in the space BV of functions 
with bounded variation, and jl0[ || for the case of L°° solutions (with possibly unbounded 
variation) of Temple class systems. 
Throughout, we shall assume: 



(Hi) the map / : £1 



i n is smooth and f2 C K n is open; 



(H2) the system (1.1) is strictly hyperbolic, i.e. the Jacobian matrix Df(u) has n real 
and distinct eigenvalues Xi(u) < ... < A„(u), u € O (with a corresponding basis of 
eigenvectors {ri(u), . . . ,r n (u)}); 

(i/3) each characteristic field ri is linearly degenerate or genuinely nonlinear in the sense of 



Lax 1 15 



(H4) there exist pe {1, n} and Co > such that 

Xi(u) < ... < X p (u) < -c < < c < X p+ i(u) < ... < A„(w), u £ Vt. (1.5) 
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By (H4), for a solution denned on the strip t > 0, x €]o, 6[, there will be n - p 
characteristics entering at the boundary point a; = a, and p characteristics entering at x = b. 
The initial-boundary value problem is thus well posed if we prescribe n — p scalar conditions 
at x = a, and p scalar conditions at x = b (see |l6f ). 

Definition 1.1 Given cp 6 i 1 ([a,6]), w £ 51, and T > 0, we sa?/ that the problem (1.1 )-(l.i ) 
is exact controllable at time T to the state v if and only if there exist measurable maps on, Pi 
such that the solution of j[l. satisfies 

u(T, •) = v, a. e. in [a, b]. 



Definition 1.2 Given (p 6 L ([a,b]), v E fi, we say that the problem (l.l)-(l.i ) is asymp- 
totic stabilizable near the state v if and only if there exist measurable maps cti , [3i such that 
the solution of \1. satisfies 

u(t, •) — > v, in ^([a, b\) as t — > +oo. 

In this note we present some recent results obtained by the authors |3| ^) concerning 
both the problem of exact controllability and of asymptotic stabilization near a constant 
state. We will first consider the case of Temple systems [E0|, for which it was obtained a 



characterization of the corresponding attainable set (L4) in terms of suitable Oleinik-type 
estimates, which is a natural extension of the results in j^, |(| concerning scalar conservation 
laws. For general nonlinear systems, one cannot expect such an exact controllability result. 
Indeed, even if all wave-fronts in the initial data exit from the interval [a, b] within finite time, 
they can generate new waves by interacting among themselves. In turn (figure 3), further 
interactions can produce a sequence of wave-fronts remaining within the interval [a, b] for 
all times t > 0. Therefore, the effect of the initial data on the solution u(T, •) may never 
be completely erased, no matter how large we choose the terminal time T. Hence, we will 
present a result concerning the asymptotic stabilization of a general system of conservation 
laws near a constant state. Finally, we discuss a counterexample to the exact controllability 
concerning a class of 2 x 2 systems for which, in general, a constant state u* cannot be 
attained, in a finite time T. 

An outline of these results established in [|[ ^) is given in the following sections. 



2 The Attainable Set for Temple Class Systems 

Our first result is concerned with the problem of exact controllability of Temple class 
systems |19L E0[, which are systems that satisfy the following additional assumption. 
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{H5) There exists a complete set of Riemann coordinates w — (wi, . . . , w n ) such that each 
level set {u ; Wi(u) = constant} is an hyperplane. 

As a consequence, all integral curves of the eigenvectors are straight lines and coincide with 
the Hugoniot curves. We shall also assume 

(Hq) as w ranges within the product set T = [w^ , w^] x • • • x [w~,w^], the corresponding 
state u remains inside the domain Q and each characteristic field is genuinely nonlinear. 

In the case of systems of this type, the well-posedness theory for the mixed problem was 
established in Q within domains of L°° functions (with possibly unbounded variation). Here, 
the boundary condition is formulated in terms of the strong L 1 trace of the solution u at the 
boundary and, in the same spirit of jl3[ ] , is based on the definition of a time-dependent set of 
admissible boundary data, that is related to the notion of Riemann problem. Moreover, for 
sake of uniqueness, it was introduced in |j, ^| a definition of entropy admissible weak solution 
to the mixed problem that includes an entropy admissibility condition of Oleinik type. 

Notice that, for such systems, wave interactions can only change the speed of wave fronts, 
without modifying their amplitudes. Therefore, the only restriction to boundary controlla- 
bility is the decay due to genuine nonlinearity. We thus consider a set of maps, defined in 
terms of the above Riemann coordinates, that satisfy certain Olcinik-type conditions on the 
decay of positive waves. 



^GL~([o,6],r) 



m(ip{y)) - Wi(tp(x)) p I for a - e - a<x <y <b, 
V-x ~ x-a I tf j e { p + i r .. )n } 

Wi(ip(y)) - Wi(ip(x)) p J for a - e - a<x <y <b, 



V~ x <>-y [if j G 

(2.6) 

The inequalities in ( |2.6| ) reflect the fact that positive waves entering through the boundary 
at a or at b decay in time. Therefore, their density is inversely proportional to their distance 
from their entrance point on the boundary. 

We can now state our first main result (see ||). 

Theorem 2.1 Let ( [_?. j[ ) be a system of Temple class, and assume that (Hi)- (Hq) are ver- 



ified. Then, letting A(T) be the attainable set defined in (in which the solution is 



understood as an "entropy admissible weak solution" ) , the following hold: 
(i) for every fixed f > 0, there exists p = p(r) > such that 

A(t)CK p , t>t; (2.7) 
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(ii) there exist T > and p' < p(T), such that 

K p ' C A(t), t>T; (2.8) 

(Hi) A.(T) is a compact subset o/L 1 ([a, b}) for each T > 0. 

The first property (i) is an immediate consequence of the definition of entropy admissible 
solution that satisfies suitable Olcinik-type estimates. 

The proof of (ii) is established in two steps. 
1) Backward Construction of Front Tracking Solutions. We take 

r > T = 4 A min = min IXA , (2.9) 

\ min ' i 

and, for a given function ip G K p , we construct a sequence of approximate solutions u v on 
the strip [0, r] x [a, b] such that 

L 1 

u u (t,-) =ip v Vv — ► V, ^ 

U u (0, •) = !!„ U„ > 93, 

with the following procedure. We partition the strip [0,t] x [a, b] in three regions. On the 
rectangle [(3/4)T, r] x [a, 6], starting from t = r, we construct backward in time the front 
tracking solution u v relying on the fact that the Oleinik estimates of the definition ([^1]) of K p 
guarantee that two rarefaction fronts of the same family never cross in Q (see figure la). The 
total number of wave- fronts in u"(t, ■) decreases as t J, (3/4)T whenever a (backward) front 
crosses the boundary points x = a, x = b. Therefore, since the maximum time taken by 



fronts of u v to cross the interval [a, b] is (b — a)/A mm , the definition (2J3) of T guarantees 
that all the (backward) fronts of u v will hit the boundaries x — a, x = b within some time 
t' G ](3/4)T, t[ . Hence, there will be some constant state u> G f2 such that ^^((3/4)^ •) = to. 
We next define u v on the rectangle [0, T/4\ x [a, b] as the restriction to [0, T/4] x [a, b] of the 



front tracking solution to the Cauchy problem for (1.1), with initial data 



u v (a+) if x < a, 
u(x) — < u v (x) if a < x < b, 
u v (b-) if x > b. 



Since u u contains only fronts originated at the points of the segment {(0,x) ; x G [a, 6]}, 
because of ( |2.9| ) these wave-fronts cross the whole interval [a, b] and exit from the boundaries 
x = a, x = b before time T/4. Hence, there will be some state lo' G VL such that ^(T/A, ■) = 
uj'. Finally, we define u v (t, ■) for t G [T/4, (3/4)T] so that u v (T/A, ■) = w', u u ((3/4)T, ■) = lu. 
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figure la 



figure lb 



2) Convergence of Front Tracking Solutions. Since the sequence of approximate solu- 
tions constructed by our algorithm satisfy the Oleinik-type estimates on the decay of positive 
waves, they have a uniformly bounded variation on each interval of the type [a+e, b— e], e > 0. 
Thus, applying Helly's Theorem, and by a diagonal procedure, we can extract a subsequence 



that converges in L 1 to a weak solution u of (1.1). We then extend a regularity property 
established by A.Bressan and P.G. LeFloch j[l| for solutions with small total variation of 
general genuinely nonlinear systems, to the case of solutions with arbitrarily large variation 
of genuinely nonlinear Temple class systems. This property guarantees that, for Temple sys- 



tems, solutions of the mixed problem (1. !)-(!. 3) are continuous outsided a countable number 
of Lipschitz curves. As an immediate consequence we deduce that the solution u admits a 
strong L 1 trace at the boundaries x = a, x = b and satisfies the corresponding boundary 
conditions. 

Concerning (iii), the compactness of A(T) is achieved with the same type of arguments 
used to establish the convergence of the approximate solutions in the proof of (ii). 



3 Asymptotic Stabilization 



We now consider a general n x n system and show that, starting with an initial data with 
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small oscillations, the system can be asymptotically steered to any constant state (see ^ 
Theorem 1]). 



Theorem 3.1 Let K be a compact, connected subset of the open domain Q C R™ and 



assume that (1.1) satisfies (Hi) — (H$). Then there exist constants Co, S, n > such that the 



following holds. For every constant state u* G K and every initial data (p : [a, b] 1— > K with 



Tot.Var.{ip} < 5, there exists an entropy weak solution u — u(t,x) of (l.l)-(l.i) such that, 
for all t > 0, 

Tot.Var.{u(t,-)} <C e~ 2Kt , \\u(t, •) - u*|| ioo < C Q e" 2 "'. (3.11) 



The idea of the proof is as follows. Call A» > a lower bound for the absolute value of all 
wave speeds and set r = (6 — a) / A* . In this way, all waves present in the solution at a given 
time t will exit through one the boundaries within time t + r. On the first interval [0, t] we 
let all waves exit, arranging the boundary values at x — a and at x = b so that no reflected 
waves ever enter the domain [a, b]. 

Therefore, the only waves present in the solution at time r are those generated by inter- 
actions, in the interior of the interval ]a, b[. They can be estimated as 



Tot.Var.{u(r, ■)} < C Tot.Var.{u(0, -)} 2 - 



The above estimate shows that the solution u(t, •) remains very close to some constant state, 
say . If ^ u* , we suitably change the boundary conditions, producing new incoming 
waves at (r, 6), (2r, a), achieving the bounds 



Tot.Var.{u(3r, •)} < C Tot.Var.{u(0, -)} 2 , ||«(3r, •) - u*\\ Loo < C Q \\u(0,-) - u*f L< 



Repeating inductively the same strategy in the time intervals [3r, 6r] , [6r, 9r] ... , we obtain 
the result. 
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4 A Counterexample to Exact Controllability 



An interesting question is whether one can reach exactly a constant state u* within a 



finite time T. By Theorem |2.1| , this is certainly the case for Temple class systems. In the 
final part of this paper, we show that this exact controllability cannot be attained in finite 
time, in general. 

Our counterexample is concerned with a class of 2 x 2 strictly hyperbolic, genuinely 
nonlinear systems with the property that the interaction of two shocks of the same family 
generates a shock in the other family (see || Theorem 2]). This is the case for the system 
(see @): 

Pt + (up) x = 0, 



Ui 



K 2 
7-1 



,7-1 



0. 



(4.12) 



with 1 < 7 < 3. Here p > and u denote the density and the velocity of a gas, respectively. 
Theorem 4.1 Consider a 2x2 system of conservation laws satisfying (Hi) and the following 
(H 7 ) there exist < A* < A* such that 

-A* < Ai(u) < -A* < < A* < A 2 (u) < A*, 



D\i ■ n > 0, 



D\ 2 • r 2 > 0, 
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ri A r 2 < 0, n A (Dri • ri) < 0, r 2 A (Dr 2 ■ r 2 ) < 0, 
where ri(u),r 2 (u) are the right eigenvectors of Df(u). 

Let (p £ BV(ja, 6];K 2 ) with small total variation and a dense set of shocks. Every entropic 
solution of (l.l)-(l.Z) has a dense set of shocks in u(t, •), for each t > 0. In particular, u(t, •) 
cannot be a constant. 



\3 



X2 



2/ 




figure 3 

Toward a proof, we note that the geometric properties of the system imply that 

a) the interaction of two shocks of the same family produces an outgoing shock of the 
other family 

b) the interaction of a shock with a rarefaction wave of the same family produces a rar- 
efaction wave in the other family. 

As in the scalar case we have (see |], Section 3]). 



Lemma 4.1 (Oleinik type estimate) Let u = u(t,x) be a solution of (1.1) with n = 2 
and satisfying (Hi) and (H-j). There exist k,S > such that, if 



Tot.Var.(u(t,-)) < 5, 



then 



Ui(t,y) - Ui(t,x) < - ■ (y - x), x<y, t > 0, i = 1,2, 



where u>x, io 2 are the Riemann coordinates associated to \l. \ ). 
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By the previous properties, a shock can never be completely canceled by interactions with 
rarefaction waves of the same family. Hence, it can only disappear by exiting from one of 
the boundaries x — a or x = b. However, if the set of shocks at time t — is everywhere 
dense, these shocks will interact among each other on a dense set of points in the domain 
[a,b] x [0, oof, and give rise to a dense set of new shocks. One can arrange so that the 
total strength of these shocks quickly approaches zero, according to Theorem 3.1, but cannot 
become exactly zero within finite time. For all details we refer to J^, Section 3]. 

Remark 4.1 The previous analysis breaks down in the case of the p— system, because in 
this case the interaction of two shocks of the same family produces a centered rarefaction 
wave of the other family. In particular, an Oleinik type estimate cannot holds. 
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